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La presente invention conceme une methode pour determiner des zones d'un milieu poreux 
10 stratifie dans lesquelles le fi'ont ou interface entre des fluides en place et des Guides de 
balayage se deplace de faqon stationnaire, c'est-i-dire sans deformation et a vitesse constante. 

La methode trouve des applications notamment quand on a a balayer par des fluides injectes 
(de Teau pai' exemple) un gisement d*hydrocarbures stratifi6. Cette determination permet 
d'obtenir ti'es rapidement les zones pouvant etre agregees et constituer des unites. 
15 hydrodynamiques lors de la mise h. Techelle d'un modele de reservoir, 

Etat de la technique : 

La simulation des 6coulements diphasiques en milieu poreux heterogene peut requ6rip5 
d'importantes ressources informatiques, en particulier lorsque le modele numerique du noilieu'- 
consider^ est fortement detailld. C'est le cas notamment en ingenierie de reservoirs, dans le" 

20 domaine p6trolier. Afin de pouvoir proceder a la simulation d'^coulements avec des moyens 
raisonnables, il faut utiliser une description reduite du reservoir. Une technique consiste k 
agr6ger les mailles du modele numerique de fa?on k obtenir un maillage plus grossier, 
compose de mailles* aux valeurs effectives moyennes. Cette technique, dite des pseudos- 
fonctions, a initialement ete propos6e pour des reservoirs stratifies de fa§on k agreger les 

25 mailles dans la direction verticale par 

- Coats, K.H., Nielsen, R,L, and Terhune, M.H. ''Simulation of three'dimensional, two- 
pJmse flow in oil and gas reservoirs" SPE 1961 , 1967 

Cette mdthode, reprise et 6tendue par de nombreux auteurs, repose sur des hypotheses 
d'equilibre vertical capillaire ou d'^quilibre gravitaire oil les forces capillaires ou gravitaires 
30 sont pr6dominantes devant les effets dus k la viscosite, et la distribution des saturations dans 
le reservoir peut 6tre connue sans recourir k une simulation sur le maillage precis (simulation 



fine). On peut en tirer facilement des valeurs moyennes suivant la direction verticale. Une 
seule couche aux valeurs moyennes peut finalement suffire k d&rire revolution des 
saturations dans le reservoir. 

Une autre methode conduisant a obtenir des valeurs moyennes est proposfe par : 

- Heam.- CX. "Simulation of stratified waterflooding ■ by pseudo-relative permeability 
curve: " Journal of Petroleum Technology, pages 805-813, July 1911. 

Elle repose sur I'hypoth^se d'dquilibre vertical du aux effets visqueux. Dans ce cas, les 
r66quilibrages visqueux dans le reservoirs sont tres xapides par rapport aux variation's de 
saturation, de telle sorte que des valeurs moyennes peuvent toujours 8tre calculees sans 
recoutir k la simulation fine, car les fluides se propagent a vitesse constants dans chaque 
couche. La mdthode de Heara, ult6rieurement Vendue pax. entre autres par : 

- Simon, A.D. and Koederitz, LJF. "M improved method for the determination of pseudo- 
relative permeability datafor stratified systems." SPE 10975, 1982, 

se base sur un calcul it6ratif conduisant k ordonner les couches du reservoir en fonction de la 
Vitesse de propagation des fluides dans celles-ci. On peut montrer que ce calcul it^ratif peut 
ne pas converger. On est alors en presence d'un front stationnaire dans au moins deux couches 
du reservoir. Les auteurs proposent alors de ne prendre qu'une seule couche aux propri6tes 
moyennes pour ces zones partjculieres. Une etude theorique de I'equilibre vertical et de ses 
implications peut 8tre trouvee chez : 

Yortsos, Y.C. "AnalyHcal studies for processes at vertical equilibrium." SPE 26022, 1992. 

Lorsque la capiUarite ou la gravite ne peuv«it 6tre considerees comme piedominantes sur les 
effets visqueux, mais que I'equilibre vertical visqueux ne peut se mettre en place, les 
methodes proposees jusqu'k pr&ent ne peuvent se passer d'une simulation fine sur tout ou 
partie du reservoir pour calculer les proprietes moyemies pouvant Stre affectees aux blocs 
grossiers. On parle alors de pseudo fonctions dynamiques. Celles-ci furent inHoduites 
notamment par : 

- Kyte. J.R. and Berry, D.W. "New pseudo fiinctions to control numerical dispersion " SPE 
5105, 1975. 
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Ces m6thodes dynamiques posent toutefois de nombreux problemes th6oriques et pratiques, 
Des revues de ces methodes et des problfemes associes peuvent Stre trouv6es dans les 
publications suivantes : 

Archer, R. ''Pseiido function generation''. Master of Science thesis, Departement of 
Petroleum Engineering of Stanford University, 1996. 

Ahmadi, A. "Utilisation des propriitis iquivalentes dans les modules de reservoir : cos 
des icoulements diphasiques incompressibles." Thfese de Doctoral, Universit6 de 
Bordeaux 1, 1992 ; 

Barker J,W. and Thibeau S. critical review of the use ofpseudo relative permeabilities 
forupscaling". SPE 35491, 1996, 

Les methodes dynamdques ont en general recours a une agregation systematique et uniforme 
des mailles. Afm de prendre en compte les influences locales des heterogeneites siir 
r6coulenient, des methodes d'agregation non-uniformes ont ete proposees notamment par : • 

- Darman, NIH. and Durlofsky, L.J. "Upscaling immiscible gas displacements: Quantitative 
Use of Fine Grid Flow Data in Grid Coarsening Schemes^' SPE 59452, 2000. 

Ces mdthodes ont Tavantage de n'agreger que des zones preferentielles et laissent en ^tat 
initial les parties du reservoir k baracteres hydrodynamiques plus complexes. EUes ne peuvent 
cependants*affranchird' une simulation fine, > 

Des methodes pemiettant d*agr6ger de fagon non unifoimes les strates d'un reservoir, sans 
avoir recours a une simulation dfitaillee, ont ete propos6es. On pent citer des m6thodes 
entiferement statiques, comme : 

- Li, D, and Beckner, B. ^'Optimal uplayering for scaleup of multimillion-cell geologic 
models'\ SPE 62927, 2000 

ou basees sur les r6sultats d'une simulation d'6coulement monophasique (qui est beaucoup 
rapide qu'une simulation d'ecoulement diphasique), comme : 

- Stem, D. and Dawson, A.G. "A technique for generating reservoir simulation grids to 
preserve geologic heterogeneity'', SPE 51942, 1999. 

De telles m6thodes, par definition, ne prennent done pas en compte les effets dynamiques et 
visqueux des icoulements diphasiques. 
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L-evo]u.on du front dans le reservoir au cours de recoulement est consid^rablement 
anfluenc^e pa,- le couplage visqueux entre le champ de pression et le champ de saturation En 
particulier, lorsque le fluide injecte est moins visqueux et par cons^uent plus mobile au 
niveau du front que le fluide en place, les instabilitds visqueuses favorisexont toujours 
I'ecoulement des fluides dans les couches les plus perm6ables. Le temps de percde au travers 
de celles-ci est beaucoup plus rapide que dans le reste du reservoir. Au contraire, si le fluide 
injects est mdins mobile, le coupfage visqueux peut ie ralentir dans les couches initialement 
plus rap,des, coxnpensant ainsi les differences de permeabiUte dues k la stratification. Un front 
stationnaire apparait alors. 

Si. dans au moins une partie du x^ervoir; 11 existe un tel front stationnaire, tout se passe 
comme sr I'ecoulement n'y voyait quHme seule couche au sens hydrodynamique. On peut done 
raisonnablement penser qu'une seule couche aux propriaes moyennes dans cette zone suffit ^ 
la moddhsation lors de la simulation num6rique de I'ecoulement. Ceci est d'importance 
capitals lorsque I'on passe du module gdologique de r^ervoir avec toutes les couches au 
15 modele de simulation, pour lequel seules les zones les plus influentes sur l'6couleme„t 
suffisent a la description. 

Deux probl^mes se posent alors. Connaissant la viscosite du fluide inject^, il faut pouvoir 
detemnner a priori I'existence et la localisation de zones ou I'ecoulement est stationnaire, sans 
recounr pour autant k une simulation complete d'ecoulement sur le module d6tailI6 n faut 
20 6galement etre en mesure avec la m6me contx^nte. de choisir la viscosity du fluide injecte de 
. telle sorte qu'un maximmn de couches aient un comportement stationnaire. 

La methode que nous proposons pennet de reconnattre les zones d'un r^servdir stratifi6 se 
comportant comme une seule couche au sens hydrodynamique sans x^couxir k une simulation 
fine et sans invoquer les hypotheses d'equihbre vertical. Hie facilite le choix des zones ^ 
25 agreger dans les reservoirs stratifies. Elle pexmet done de prendre en compte les effets 
dynam^ques et visqueux tout en permettant une determination des couches gxossi^res tr^s 
rapide par rapport aux anciennes solutions. 

La methode selon I'liivention 

U methode selon I'mvention pennet de d^tenniner dans un nulieu s,ratifi6 dont les profriitM 
30- phys,ques son, oonnues ou esMmaes, au moins une zone oD une interface entte un fluide en 
plaoe dans le milieu et un fluide de balayage. de viscosity et de„a,t& difK^tes connues 
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inject&s dans le milieu, se d6place de fagon stationnaire, dans le but de simplifier la 
construction d'un modele de simulation des ecoulements dans le milieu. EUe comporte 
essentiellement les 6tapes suivantes : 

- on considdre une forme d'interface a priori en supposant que les d6placements de fluides 
5 en tous ses points sont stationnaires ; 

- on d6tennine le champ de pression de part et d'autre de Vinterface a priori ; 

- on modifie it6rativement la forme de Tinterface jusqu'a obtenir que les pressions de part et 
d*autre d'au moins une partie de Tinterface s*6galisent en tous points de cette partie ; et 

- on attribue des proprietes hydrodynamiques moyennes unifoimement & chaque zone du 
10 milieu delimitee par chaque partie dlnterface, quand la dite egalisation est atteinte. 

Suivant un mode de mise en oeuvre, quand on ne pent obtenir une egalisation des pressions de 
pait et d'autre de Tinterface tout le long de celle-ci, on segmente interface en plusieurs parties 
et on modifie iterativement et separement la forme de ces differentes parties, jusqu'k obtenir 
une Egalisation des pressions de part et d'autre de chacune d'elles, Tetendue de chaque partie 
15 d'interface, quand la dite egalisation est atteinte, delimitant une zone favorable a laquelle on 
attribue uniformfiment des propri6t6s hydrodynamiques moyennes. 

Suivant un exemple de mise en oeuvre, on determine la forme d'au moins une zone du milieu 
delimitee par une interface k deplacement stationnaire, qui correspond k des valeurs 
differentes de la viscosit6 des fluides de balayage, et on choisit la viscosite pour laquelle on 
20 optimise les deplacements stationnaires dans le dit milieu. 

Presentation succincte des figures 

Les caractdristiques et avantages de la m6thode selon Tinvention, apparaJtront plus clairement 
k la lecture de la description ci-aprfes d*un exemple non limitatif de naise en c&uvre, en se 
referant aux dessins annexes oii : 

25 - la figure 1 montre une boucle de determination iterative de Texistence ef de la forme dun 
front stationnaire pour un dcoulement diphasique de rapport de viscosites connu dans un 
milieu stratifie donnd ; 



- la figure 2 niontre une boucle de construction itemtive de revolution- des couches 
hydrodynamiques en fonction du rapport des viscosity, pour des 6coulements diphasiques 
dans un milieu stratifie donn6 ; 

- la figure 3 montre un exemple de i^servoir stratifie incline par rapport k la verticale, Cette 
5 figure a pour but d'expliciter les notations employees dans la description de la m6thode ; 

' ^ ^^""^'^ ^'^ perm^abilites absolues en fonction de la profondeur " 

employees dans I'exemple synth^tique de reservoir iUustrant la m6thode de fafon non- 
limitative ; 

- la figu,^ 5 montt^ revolution de I'interface au cours du processus iteratif de recherche de 
10 la fonne d'un front stationnaire, sur I'exemple decrit par la figure 4 ; 

- la figure 6 montre revolution de la solution stationnaire en fonction du rapport des 
viscosites employees, sur I'exemple d6crit par la figure 4 ; et 

- la figure 7 montre revolution du nombre de couches hydrodynamiques en fonction du 
rapport des viscosites employees, sur I'exemple decrit par la figure 4. 

15 Descripaon detainee de la methode 

On consid^re un balayage de Guides de viscosites connues dans un milieu stratifie. Supposons 
qu'un front stationnaire de forme connue se developpe dans le miUeu consider^ Ceci impose 
une condition de flux au niveau de I'interface. Le flux est egalement connu sur les bords du 
reservoir (o5 n est en general nul). en entire et en sortie (ot il est impose par les conditions 
20 d'mjection et de recuperation). Comme la pression est une variable continue, le saut de 
piession de part et d'autre de I'interface est nul en tout point de ceUe-ci. 

On cherche k present I'existence et la forme d'un front stationnaix^ dans le milieu. On procMe 
de fagon iterative en partant d'une forme de front a priori, en supposant que celle-ci 
coxrespond k la solution stationnaire. Comme le front est suppose stationnaire, on connait le 
25 flux k I'interface et I'on a done un probl^nie de Neumann parfaitement defmi k resoudre pour 
chacune des deux zones que separe I'interface. champ de pression pent done etre resolu 
.ndependamment dans les deux regions. On evalue alors le saut de pression de part et d'autre 
de I'interface, en tout point de celle-ci. Si le saut de pression est nul en tout point, on poss^de 
alors la forme coirecte du front conespondant k un ecoulement stationnaire. Si le saut n'est 
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pas nul partout, la condition de continuite des pressions n'est pas respect6e et la forme de front 
ne correspond pas k une solution stationnaire. On doit alors deformer Tinterface et 
recommencer le prdcessus jusqu'a la convergence. Deux cas apparaissent alors : 

- Si la convergence est atteinte, on possdde alors la forme du front correspondant a une 
5 solution stationnaire de l'6coulenient. Toutes les couches sont travers6es par le fluide a la 

meme vitesse, et Tensemble du systeme se comporte comme une seule unit6 
hydrodynanoique. On peut done ne considerer qu'une seule couche aux proprietes 
moyennes lors de la naise a Techelle, 

- Si le processus it6ratif ne converge pas, c'est qu'un front stationnaire ne peut pas exister 
10 sur Tensemble du milieu. Lors de la deformation du front, il apparait alors en general deux 

zones s'eloignant perpetuellement Tune de Tautre. On decoupe alors le milieu en deux 
parties correspondants aux zones pr6cedemment identifiees. Puis le processus de 
recherche d'une solution stationnaire reprend dans chacune de ces deux parties. En 
prccedant ainsi de maniere iterative, on finira par obtenir toutes les zones dans lesquelles 
15 des fronts stationnaires peuvent se d6velopper, et qui peuvent Stre consid6rees comme des 

unit6s hydrodynamiques. 

La m6thode peut aussi etre appliquee de fagon iterative eni faisant varier le rapport des 
viscosit6s des fluides consideres. On part d'un fluide injects beaucoup plus visqueux que celui 
en place. On doit alors trouver rapidement une solution stationnaire sm- toute r6paisseur du 
20 milieu, Lorsque Ton dinainue le rapport des viscosites entire le fluide inject6 et le fluide en 
place, la forme du front stationnaire, si celui-ci existe, doit 8tre modifi^e. On reprend alors la 
methode dficrite ci-dessus pour trouver la nouvelle fomie du front, en utilisant connme forme 
de front initiale du processus iteratif, la forme finale, obtenue pour le pr6c6dent rapport des 
mobilit6s. Ceci permet de converger beaucoup plus rapidement vers la solution stationnaire. 

25 On peut ainsi proc6der de fagon iterative sur le rapport des viscosites- Au fur et k mesure que 
Ton d6crott ce rapport, le couplage visqueux ne sera plus suffisant pour compenser les 
contrastes de perm6abilite et une solution stationnaire ne pourra plus Stre trouv^e sur toute 
Tepaisseur du milieu. On d6coupe alors le milieux en deux portions suivant la m6thode decrite 
ci-dessus en cas de non-convergence, et la recherche se poursuit dans les deux portions de 

30 milieu. En procedant ainsi de fagon it6rative sur le rapport des viscosit6s, on obtient le 
nombre de couches hydrodynamiques du noilieu en fonction des viscosites consid6rees. Ce 
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r6su]tat pourrait permettre par exemple de connaJtre la viscosite du fluide d'injection qui 
permet d*optimiser la recup6ration. 

On consid&re une section verUcale d'un reservoir stratifie selon une direction x inclinde par 
rapport ^ rhcrizontale (figure 3). Les valeurs des permeability absolues sont connues dans 
5 toutes les strates. L'epaisseur du reservoir est selon la direction z. A rextr6mit6 amont {x = 0) 

°" ^ debit constant j2o.impos6, 4an fluide incompressible de-viscosif^ et de densite 

Pu qui pousse de fagon immiscible le fluide initiaJement pi^ent dans le i^ervoir 
mcompressible et de viscosity et de density p.. On suppose qu'il n'y a pas de pression' 
capillaire entre les deux fluides. La r&up^rBtion se fait h rextr6mit6 aval (;c = L) Le flux est 

10 nul a trovers le toit (. = 0) et le mur du i^ervoir .= A). On suppose egalement que les 
courbes de penneabilit6s relatives sont unifonnes dans le resewoir et que chacune de ces deux 
zones est i saturation constante en fluides consideres. La zone non encore balayee est done ^ 
saturation itrdductible en fluide injecte et la zone balayee ^ saturation maximum ou 
saturation de front Sf selon les hypoth^es de I'utilisateur, Le front de saturation est 

15 arbatrairement place au centre du reservoir, avec une forme px^suppos6e qui conespond a 
notre a pnori sur la foime du front stationnaii-e. Cette forme sera par la suite modifi6e jusqu'S 
trouver la "vraie" solution stationnaire, si celle^^i existe. On doit done :€soudie un probl^me 
aux frontidres libres. 

Les conditions aux limites et la condition de front stationnaire pemiettent de definir un 
20 probl&me de Neumann dans les deux zones que separe le front : 

- A I'extremit^ amont : la condition d'injection impose un flux constant sur l'epaisseur du 
r6servoir en entirde. Si le miUeu est suppose suffisamment long et le profil de saturations 
deja bien ddvelopp6. on pent supposer que le gradient de pression en entrde du syst^me est 
identique dans toutes l^s nr.n/->,*.c • 



identique dans toutes les couches : 
25 dp 

ax 



- ^>^o ^ ^ „ 
p^ + Ag-x 



dp 

dz 



jr»0 



cap «t la pKssion au poin, consid^rf e. if. est la moyen„e des permlabilit& effecUves ™ 
«.Wc sar «>ute r^paisseur du „»li«,. « fi. es, Ic d^bi. d'tojecdon. Lea p=™&bili«a 
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effectives prennent en compte Teffet des penneabilites relatives, en fonction des saturations 
choisies par Tutilisateur dans les zones amont et aval, 

- A Vextremite aval : la condition d'incompressibilite impose que le flux sortant du systeme 
soit le meme que le flux entrant. On obtient done une condition similaire a la tranche de 
. recuperation : 



dx 



dp 

dz 



ou K2 est la moyenne des penn^bilit6s effectives en sortie sur toute r6paisseur du milieu. 

- de part de d'autre de Vinterface : Comme le front est suppos6 statiotinaire, il se d6place 
horizontalement sans d6formation k la vitesse de filtration U„ , d6finie telle que : . 



X 



15 



20 



Oh X est le vecteur unitaire dans la direction horizontale. Cette condition impose le flux et 
done le gradient de pression de part et d'autie de I'interf ace : 



JEl 



Ml M2 

oH nest le vecteur xinitaire orthogonal h. I'interface au point consid^re, Vpj est le gradient de 
pression h I'interface du c6t6 amont et Vpj est le gradient de pression h. I'interface du c6t6 
aval. 

- Equation en pression : L'equation en pression est donnee classiquement de part et d'autre 
de I'interface. Pour la zone aval : 



r 



Ml 



(Ag-Vp) 



= 0 
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Pour la zone amont 




Les equations en pression et les conditions aux Hmites pennettent de d6fmir pour chaque zone 
de part et d'autre du front un probldrae de Neumann. Le champ de pression dans chaque zone 
pent done Stre dalcule num6iiquement h une constante pres pour chaque zone. On limine ce 
problfeme en posant la condition suppldmentaire suivante, pennettant d'eUminer les constantes 
en recalant les pressions le long du front : 

otpj et p2 sont les pressions 6valu6e au front respectivement pour les zones amont et aval, et 
C est d6finie arbitrairement par I'utilisateur et peut atre nuUe. 

Une fois la pression p calcuMe dans chaque region et recal^e le long de I'interface. le saut de 
pression Sp^ (z)k travers le front en tout point de ce dernier peut Stre calculi : 

^Pf (2) = Pa Cz) -Pi (z) - Ayo(g.x + g.z) 

Si le saut de pression est partout nul le long du front, ou du moins en de?^ de la valeur 
critique ddfinie par I'utilisateur, on dispose de la forme du front coirespondant k un 
d^placement stationnaire. Sinon, il faut dgfoimer interface de fa^on k faiie diminuer les sauts 
de pression jusqu'k obtenir la forme stationnaire. Ceci revient alors h i^soudre an problSme 
d'optimisation, pour lequel de nombreuses m6thodes sont applicables. 

Ce processus de minimisation peut ne pas converger. Dans ce cas, il n'existe pas de front 
stationnaire sur toute I'epaisseur du milieu. La m6thode decrite doit done 6tm reprise sur des 
parties plus petites du reservoir. On peut ainsi tirouver des fronts stationnaires dans plusieuis 
zones du r6servoir, qui se deplacent k des vitesses diff^rentes. 

Exemple 

On d€crit ici un exemple synth^tique illustirant la m^thode de fagon non-limitative. 

On considdre un reservoir stcatifi^ avec quati-e strates horizontales de permeability 
diff^rentes, indiqudes sur la figure 4. Pour simplifier au maximum, on n^gligera les effbts de ' 



la gravite et on considerera que les perm6abilites relatives des stxates sont telles que la zone 
non balayee est satur6e en huile seule et qu'en amont de Tinterface, seule Teau est presente k 
saturation 5=7. La porosite est constante dans le reservoir. Ueau est inject6e sur la face amont 
k flux constant connu. L'huile est i viscosit6 de L 

5 On commence par chercher Texistence et la forme d'un front stationnaire dans le milieu 
lorsque la viscosit6 de I'eau est de 12, Comme on ne dispose pas encore d'a priori, on part 
d'une interface plane, que Ton deforme de fagon iterative jusqu'^ atteindre la convergence vers 
une solution stationnaire (figure 5). 

Urie fois la solution obtenue pour cette viscosite, on peut rechercher la solution stationnaire 
10 pour une viscosit6 de Teau inferieure, par exemple une viscosity de 10. La solution 
stationnaire obtenue precedemment peu servir d'a priori en entree du processus iteratif de 
deformation du front. On converge alors vers une nouvelle solution. En proc6dant ainsi 
it6rativement sur difffirentes valeurs de la viscosite de Teau, on obtient les formes du front 
stationnaire associ6es k ces valeurs (figure 6). Corome le rapport des mobility entre'le fluide 
15 ne place et le fluide inject^ diminue chaque fois que Ton diminue la viscosite de Teau, le front 
s'etale au fur et k mesure, car les effets visqueux ont de plus en plus de mal k compenser les 
differences de perm6abilit6 d'une strate a Tautre. Lorsque la viscosite de I'eau diminue en degk 
de la valeur critique 5, on s'aper$oit que le processus de recherche d'un front stationnaire ne 
converge plus. Les positions du front dans les deuxieme et troisieme strates s'eloignent 
20 perpetuellement Tun de Tautre. On considere alors qu'un front stationnaire ne peut se former 
sur toute I'^paisseur du milieu. Ce dernier doit Stre decoupie en deux parties, la s6paration se 
situant entre la deuxifeme et la troisi&me couche. On poursuit ensuite le processus de recherche 
de front stationnaire dans chacune des parties pr6c6demment identifies. 

En iterant ainsi sur la valeur de la viscosit6 du fluide injecte, on peut construire un diagranmie 
25 representant la leponse du milieu du point de vue hydrodynamique en fonction du rapport des 
mobilit6s M entre le fluide en place et le fluide injecte (figure 7). On voit sur ce diagramme 
que lorsque le rapport des mobilit6s est sup6rieur k 5, le milieu se comporte conrnie une seule 
couche hydrodynamique, puisqu'un front stationnaire existe sur toute son 6paisseur. Lorsque 
M passe entre 4,5 et 5 deux couches hydrodynamiques apparaissent : la premiere est formee 
30 des deux premieres couches "geologiques", et la seconde des deux autres couches 
"g6ologiques". Puis les deux premiSres couches geologiques se s6parent au sens 
hydrodynamique en degk d'un rapport M de 4,5. On a alprs trois couches au sens 
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hydrodynamique dans le milieu. Les deux demi^res couches g^ologiques se s^parent en degk 
de M=3. On a alors quatre couches hydrodynamiques, qui correspondent aiix quatre couches 
g^ologiques prdsentes dans le milieu. 
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REVENDICATIONS 

1) M6thode pour ddterminer dans un milieu stratifie dont les proprietes physiques sont 
connues ou estimees, au moins une zone oil une interface entre un fluide en place dans le 
milieu et un fluide de balayage, de viscosites et densit6s diff6rentes connues, injectes dans le 

5 milieu, se deplace de fafon stationnaire, dans le but de simplifier la construction d*un module 
de simulation des 6coulements dans le milieu, caractfirisee en ce que : 

on considere une forme d'interface a priori en supposant que les deplacements de fluides 
en tons ses points sont stationnaires ; 

- on determine le champ de pression de part et d'autre de Tinterface a priori ; et 

10 - on modifie it6rativement la forme de I'interf ace jusqu'a obtenir que les pressions de part et 
d'autre d*au moins une parde de I'interface s'egalisent en tous points de cette partie ; et 

on attribue des proprietes hydrodynanoiques moyennes uniformement a chaque zone du 

milieu delimitee par chaque partie d*interface, quand la dite 6galisation est atteinte. 

. ' ■ * 

2) Methode selon la revendication 1, caract6risee en ce que, a defaut d'obtenir une 6galisation 
15 des pressions de part et d'autre de Tinterface tout le long de celle-ci, on segmente I'interface en 

plusieurs parties et on modifie iterativement et s6parement la foraie de ces differentes parties, 
jusqu'a obtenir une egalisation des pressions de part et d'autre de chacune d'elles, Tetendue de 
chaque paitie d'interface, quand la dite egalisation est atteinte, delimitant une zone favorable k 
laquelle on attribue uniforaiement des propri6tds hydrodynamiques moyennes. 

20 3) Methode selon la revendication 1 ou 2, caract^risee en ce que Ton determine la forme d'au 
.moins une zone du milieu delimitee par une interface k deplacement stationnaire, qui 
correspond k des valeurs diff6rentes de la viscosite des fluides de balayage, et on choisit la 
viscosite pour laquelle on optimise les deplacements stationnaires dans le dit milieu. 
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